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Introduction 



The (two-sided) fractional Brownian motion t ^ Bt, t G M (fBm for short) 
with Hurst exponent a, a G (0, 1), defined as the centered Gaussian process 
with covariance 

E[BsBt] = ^(|sp" + -\t- s|2"), (0.1) 

is a natural generalization in the class of Gaussian processes of the usual 
Brownian motion (which is the case a = ^), in the sense that it exhibits 
two fundamental properties shared with Brownian motion, namely, it has 
stationary increments, viz. E[{Bt-Bs){Bu-By)] =E[{Bt+a-Bs+a)iBu+a- 
By^a)] for every a,s,t,u,v E M, and it is self-similar, viz. 

VA > 0, {Bxt,t G M) ^ (X^Bt, t E M). (0.2) 

One may also define a d-dimensional vector Gaussian process (called: d- 
dimensional fractional Brownian motion) by setting Bt = {Bt{l), . . . , Bt{d)) 
where {Bt{i),t G IK)i=i,...,d are d independent (scalar) fractional Brownian 
motions. 

Its theoretical interest lies in particular in the fact that it is (up to 
normalization) the only Gaussian process satisfying these two properties. 

A standard application of Kolmogorov's theorem shows that fBm has 
a version with a~-Holder continuous (i.e. K-Holder continuous for every 
K < a) paths. In particular, fBm with small Hurst parameter a is a natural, 
simple model for continuous but very irregular processes. 

There has been a widespread interest during the past ten years in con- 
structing a stochastic integration theory with respect to fBm and solv- 
ing stochastic differential equations driven by fBm, see for instance [15\ 
El [H [22| I23j . The multi-dimensional case is very different from the one- 
dimensional case. When one tries to integrate for instance a stochastic 
differential equation driven by a two-dimensional fBm B = (i?(l), -6(2)) by 
using any kind of Picard iteration scheme, one encounters very soon the 
problem of defining the Levy area of B which is the antisymmetric part of 
Ats := £dBt,{l) dBt^{2). This is the simplest occurrence of iterated 

integrals B^^{ii, . . . ,ik) := jl dBt^{ii) . . . Jl'"'^ dBt^{ik) , ii,...,ik < d for 
d-dimensional fBm B = {B{1), . . . , B{d)) which lie at the heart of the rough 
path theory due to T. Lyons, see \17\ I18j. An alternative construction has 
been given in [lOj under the name of 'algebraic rough path theory', which 
we now propose to describe briefly. 
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Assume Tt = (rt(l), . . . ,Tt{d)) is some non-smooth d-dimensional path 
which is a-Holder continuous. Integrals such as f fi(rt)drt{l) + ... + 
fd{Tt)dTt{d) do not make sense a priori because T is not difFerentiable 
(Young's integral \16^ works for a > ^ but not beyond). In order to define the 
integration of a differential form along F, it is enough to define a truncated 
multiplicative functional or geometric rough path (F^, . . . , rL^/"J ) lying above 
r, [l/aj=entire part of 1/a, where F^^ = {6r)ts ■= — is the increment 
of F between s and t, and each F*^ = (F'^(ii, . . . ,ik))i<ii,...,ik<d: A; > 2 is a 
substitute for the iterated integrals /J(iFti(ii) J^^ drt2{i2) ■ ■ ■ fs''~^ dTt^.{ik) 
with the following three properties: 

(i) {Holder continuity) each component ofT^ is ka^ -Holder continuous, 
that is to say, kn-Holder for every k < a; 

(a) (multiplicative or Chen property) letting ST^^^ := Fj^ — F^^ — T^^, one 
requires 

dTl,{ii,...,ik)= Yl r,^^(n,---,%)rSlfe+i,---,4)- M 

ki+k2=k 

(Hi) (geometric or shuffle property) 

Fj^^ (ii, . . . , i„j)Fjj^(ji, . . . , jnj) = Fj^^'^ ^(/ci, . . . , fc^j+nj) 

fcesh(i,i) 

(0.4) 

where Sh(i,ji) is the subset of permutations of ii, . . . ,ini, ji, ■ ■ ■ , jn2 
which do not change the orderings of {ii, . . . , im) and (ji, . . . ,jn2)- 

Then there is a standard procedure which allows to define out of these 
data iterated integrals of any order and to solve differential equations driven 
by F. 

The multiplicativity property (j0.3p and the geometric property (|0.4p are 
satisfied by smooth paths, as can be checked by direct computation. So the 
most natural way to construct such a multiplicative functional is to start 
from some smooth approximation F^, e ^ of F such that each iterated 
integral T^f {ii, . . . k < [l/a\ converges in the /cK-Holder norm for 
every k < a. 

This general scheme has been applied to fBm in a paper by L. Coutin and 
Z. Qian [5] and later in a paper by the author [26], using different schemes 
of approximation of S by a family of Gaussian processes (living on the 
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same probability space) with e — > 0. In both cases, the variance of the Levy 
area has been proved to diverge in the hmit e ^ when a < 1/4. 

Let us explain briefly our construction for the second-order iterated inte- 
gral B|^(ii,i2) (by abuse of language, we shall call this object a Levy area, 
although the Levy area is usually defined as the corresponding antisym- 
metrized quantity). 

Let a G (0, ^) (or even a G (0, 1/4)), and consider the natural iterated 
integral {ii,i2) for some family of approximations B"^ of B. Assume 
Z'^(ii,i2) is some a.s. 2a~-Holder random function (living on the same 
probability space as -B^), antisymmetric in {11,12)- Then 

7^B2f(^l,^2) :=B,Y(ii,i2)+<5Zf,(n,i2) = B2f(ii,i2)+(Zf(n,i2)-Z^(ii,i2)) 

(0.5) 

satisfies properties (i), (ii) and (iii). The multiplicativity property is pre- 
served because 5 acting on an increment vanishes }10], 

6{6Z'{h,i2))tus = (Z?(il, i2) - Z^(ii, i2)) - (Z?(il, i2) - Zi{h,i2)) 

- {Zl{n,i2)-mh,i2)) = 0, (0.6) 

and the geometric property is also preserved because Z^(ii, 22) + Z^(z2, ii) = 
0. 

The function may be seen as a counterterm. Now (see section 1) Z^ 
may be chosen so as to make the regularized Levy area TZBff {ii,i2) converge 
for every k < a and T > in L'^{Q,C'^'^{[—T,T])) to a finite Levy area 
for B, where C'^'^{[-T,T]) is the Banach space of 2ft;-Holder 1-increments 
on [—T,T] in the sense of Gubinelli [10], equipped with the Holder norm 
||/||2k,t = supj,^^g[_j^ jtZ^pT- More precisely, Z^ may be chosen in the 
second chaos of B. One may prove bounds of the type 

E7^B,Y(^l,i2)^ < C|t-s|^°, E(7^B^2f (ii,i2) -^B,2f (ii,^2))^ < Ck-r/P" 

(0.7) 

from which (see |27) . Proposition 1.5) the convergence in L'^{Q, C'^'^{[—T, T])) 
may be proved by using standard arguments, in particular the Garsia- 
Rodemich-Rumsey lemma [8]. 

The first section concerns the construction of the counterterm Z^ and 
the regularization of the Levy area. The main tool is Fourier transform. 
Splitting iterated integrals into increment/boundary terms (see below) and 
reordering Fourier components by Fubini's theorem in such a way that in- 
nermost integrals bear highest Fourier frequencies ensures the proper Holder 
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regularity for every boundary term, hence the name of Fourier normal order- 
ing. Then increment terms are regularized by introducing an adequate cut 
in the Fourier domain of integration. The multiplicative rule is preserved 
by doing this (see eq. (j0.6p ). 

The sketchy generalization to higher-order iterated integrals in section 
2 involves tree combinatorics linked with the Hopf algebra structure of the 
algebra of decorated rooted trees, as defined by A. Connes and D. Kreimer 
in a series of papers [21 [3l Hj devoted to a mathematical axiomatization 
of the BPHZ (Bogolioubov and coauthors) algorithm of renormalization of 
Feynman graphs in quantum field theory [12]. Our algorithm has nothing 
to do with BPHZ though, since it is based on a somewhat arbitrary but 
convenient regularization of skeleton integrals (higher-order generalizations 
of the increment term of the Levy area) , which are tree iterated integrals of a 
particular type. The proof of the multiplicative/Chen and geometric/shuffle 
property in [28] relies on Hopf algebra computations though. A BPHZ or 
dimensional renormalization scheme could be applied to skeleton integrals 
instead of a blunt regularization (work in progress) . This algebraic approach 
has proved to be useful in a variety of contexts ranging from numerical 
analysis to quantum chromodynamics or the study of polylogarithms or 
multi-zeta functions, see for instance [141 1211 [29] . 

1 Definition of a Levy area 

Recall that fractional Brownian motion B may be defined via the harmo- 
nizable representation [21] 

Bt = co. [ w(dO (1-1) 

Jm K 

where {W^,(, > 0) is a complex Brownian motion extended to M by setting 

w.^ = -w^ (e > 0), and c„ = iyrz=^. 

We shall use the following approximation of -B by a family of centered 
Gaussian processes {B^,e > 0) living in the first chaos of B. 

Definition 1.1 (approximation B'^) Let, for e > 0, 

Bf = cj e-^l«l|e|^-"^^^ W{di). (1.2) 
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The process B'^ is easily seen to have a.s. smooth paths. The infinites- 
imal covariance E(i?^)^(i?^)J may be computed explicitly using the Fourier 
transform [6] 

V27r Jr 2cos7rar(-2a) 

(1.3) 

where k's' {s - t) := - t) + 2r/)2"-2. By taking the real part 

of these expressions, one finds that B'^ has the same law as the analytic 
approximation of B defined in [26] , namely, B^ = Tt+ie + ^t-ie = 2Re Tt+ie , 
where F is the analytic fractional Brownian motion (see also |25j). 

This second-order case is too simple to capture the combinatorial features 
of Fourier normal ordering. On the other hand, the proof of convergence for 
the regularized Levy area (after substraction of the counterterm) is short, 
and the proof for tree integrals in the general case may be considered as a 
generalization. 

Fourier normal ordering is the combination of (i) a number of equivalent 
splittings of the iterated integrals related by Fubini's theorem; (ii) a split- 
ting of the Fourier domain of integration into a number of disjoint Fourier 
domains; (iii) an appropriate choice of splitting of the iterated integrals on 
each Fourier domain. 

Let us first write down the two equivalent splittings in the case of the 
Levy area. The components -B(l), -6(2) are assumed to be constructed from 
i.i.d. Brownian motions VF(1), H^(2) via the above harmonizable repre- 
sentation. Note that Bj^^(j,j) = ^{Bf{j) — i?|(j))^, j = 1,2 requires no 
regularization, hence the only problematic second-order iterated integrals 
are the mixed integrals B^'^(i, j), i ^ j, say, i = 1, j = 2. 

Lemma 1.2 Lei (1, 2) = di?^^(l) di?^^(2). Then: 



/+00 r+co 
dw^,{i) / ^iT^6(2)e-^^'«^'+'^^'^iei6i^-"/<s(ei,e2) 
-oo J — oo 

(1.4) 

where ^ 

/ts(ei,6):= / e™i«idni P e'^^^^^^- (1-5) 

J s J s 

The proof is straightforward. 

There are two apparently equivalent ways of splitting the integral Its{£,i-,S,2) 
into an increment term, {5J)ts = Jt — Js, and a boundary term denoted 
by the symbol d: 
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(i) either writing 1^,(6,6) as (G+ - + where 
(provided + 6 / 0) 

(1.6) 

(ii) or (using Fubini's theorem) 

rt ft 



' «2 



/is (6, 6) = / e^"2«2d7X2 / 

J S J U2 

= {g; -G;){C2,^i) + irsiC2,Ci){d), (1.7) 

where 



(li 



In either case, the inner integral J^e^^^du', (ii,x,^) = (ni,s,^2) or 

r 



{u2,t, ^i) has been formally decomposed as J" e^'^'^du' — e^'^'^du' = 



which introduces an apparent infra-red divergence, since the single 
terms diverge when C — > 0. Boundary terms come from the 

contribution of ^j^, x = s or t, and are not increments. 

The idea of Fourier normal ordering is that innermost integrals should 
bear highest Fourier frequencies in order to get correct Holder estimates 
separately for the increment and the boundary term. Namely, using for in- 
stance the decomposition eq. (11.60 for arbitrary values of ^1,^2 yields a 
boundary term —B^{2){B^{1) — B^{1)) which is obviously only a~-Holder, 
and not 2a^-H61der. 

We shall say that a function of two arguments (t, s) — > Xj^ (e > 0) in 
the second chaos of B is uniformly 2a~ -Holder in e if 

(i) E(Xf^)^ < C|t — for a constant C which is independent of e, and 
furthermore 

(ii) one has the following rate of convergence E (Xf^ — Xj^)^ < C|e — 
when e,r] ^ 0. 



^Formally J" e^^^du = 4^ — j^^^ e^^^du depending on the sign of ^. 
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This implies (by the arguments given in the Introduction) that X^^ con- 
verges in L2(17;C2'^([-T,r])) for every T > and k < a. 

Lemma 1.3 (i) (boundary term^ For every a E (0, |), 
B?f+(l,2)(a) -.= 01 



f [ dW5,(l)dt^5,(2)e-^(l«^l+l«^l)|ei6l^-"4(6,6)(5) 
•''I6I<I6I 

(1.9) 



and 



/ / (iTyg,(i)(iTy5,(2)e-^(i«^i+i«^i)|66l^-"/i;(6,ei)(5) 

•''I6I<I6I 

(1.10) 

are 2a -Holder uniformly in e. 
(a) (increment termj For every a £ (1/4,1/2), the functions 

<5G?f+(l,2):=c2 / / e-(l«^l+l«^l)|6e2|^-"(G'+(ei,6)-G+(6,6)) 

(1.11) 

and 

6Gt~il,2):=cl [ [ e-(l«^l+l61)|^^^2||-(G-(e2,ei)-G;(e2,6)) 

(1.12) 

satisfy the first estimates (i) El^G^^f '=^(1, 2)p < C\t - s\^'^. 
Remarks. 

1. Only the increment 5G^'^'^(1,2) makes sense: ^^'"^'^(1,2) defined as 
the integral of G^(^i,^2) is infra-red divergent (see also remark after 
Lemma 11.41 below) . 

2. Evenifa > 1/4, only the sum of the two increment terms 6 G^'f'~^ (1,2) -\- 
6Gff'~ (1,2) satisfies the above rate of convergence estimate (ii). This 
is due to the spurious singularity on the diagonal = —^2, not to an 
ultra-violet divergence when — *■ 00. We skip the proof which 
is not needed. 

Proof. Let us first remark that the symmetry (VF(1) ^ W{2), s ^ t) ex- 
change_s ^+(6, 6)-G+(6, 6) with G; {^2,Ci)-G~ {^2,^1), and /+(6, 6)(9) 
with I^^ {C2,(,i)id)- Hence it is enough to prove Holderianity for B^''^'+(l, 2){d) 
and 5G2'^'+(1,2). 
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We shall use a number of times the following elementary lemma, inspired 
by arguments of J. -P. Kahane concerning the regularity of random Fourier 
series [131: 



Lemma 1.4 (i) Let F{u) = f^dW^aiOe'''^ , where |a(OP < C\^\-^-^^ 
for some < /5 < 1; then, for every ui,U2 £ M, 

E|F(ui) - F(u2)|^ < C'\ui - U2p^. (1.13) 

(ii) Let F{e) = J^dWi.a{Oe~''\^\ {e > 0), where |a(OP < C\C\~^-^^ for 
some < (3 < 1: then, for every ei,e2 G M+, 

E\F{e,)-F{e2)\^<C'\e,-e2\'"'. (1.14) 

Proof. Bound |e™i€ - e^^^q _ ^^^H^l for |^| < and by 2 

otherwise, and similarly for le""^!'^' — e~^2'^'|. Note the variance integral is 
infra-red convergent near ^ = 0. □ 
Remark: Unless |a(OP is Lj^^ near ^ = 0, only the increments F{ui) — 
F{u2), F{ei) — F{e2) are well-defined. 

(i) Apply Lemma 0(1) to K,(n) = ^ dW^^ (l)a(Ci)e'<i with 



•^l6l>ICi| 

since Var a(^i) < Clc^il"^"^", one gets the uniform Holderianity es- 
timates (i) E(Bj2f'+(l,2)(a))2 < C|t - for Bff'^{l,2){d). As 
for the rate of convergence (ii) (see above Lemma II. 3p , one rewrites 

e-£(iai+iai) _ e-^(l€il+l«2|) as (e"^l^il - e"''l^il)e"^l?2| (,-v\ii\^g-e\i2\ _ 

g-r?[C2|^^ The first term may be bounded as Fg by applying Lemma 
11.41 (ii). For the second term, Lemma 11.41 (ii) should be applied to 
Fs{e) = (iW^2(2)e~^l^2|gisC2|^2|~2~", which yields an exponent 

2a instead of 4a. 

(ii) Apply Lemma O (i) to F{u) = dM/ga(Oe'"^ with (setting ^ = 

ei + 6) 



a(e) = -/ dWi:,\^-C2\"''—, e-^^l*2|+i4~42U_ (i^ig) 

1? J\i2\>\i-i2\ 1« 



Setting e directly to yields (assuming for instance > 0) 

1 r+oo 

Var a{0 < / |6r'-'°lC - ^2\'-"'d^2, (1-17) 
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which converges if and only if a > 1/4, in which case 

r+oo 

Var a(e) < |er^~^° / |n|-i-2a|i _ u\^--2»du = C7|Cr^-^". (1.18) 

□ 

All together one has proved (up to the proof for the rate of convergence 
for that 

= (B,Y'+(l,2)(a) +(1,2)) + (B?f-(l,2)(a) +5G?f-(l,2)) 

(1.19) 

is 2a -Holder uniformly in e, provided a > 1/4. 
So what should one do when a < 1/4 ? 

Definition 1.5 (cut Fourier domain) Let, for some constant Creg G (0,1), 

Keg -={{^1,^2) I 161 < 1^21,16+^21 >ae9le2|}. (1.20) 

The condition |6 +61 > Creg|6l excludes a conical region along the 
singular line 6 = ~6- 

Definition 1.6 (regularized Levy area TZ^ff ) Let 

7^B?f :=(B,2f+(l,2)(a)+<57^G,2f'+(l,2))+(B,Y'-(l,2)(a)+^7^G,Y'-(l,2) 

(1.21) 

where the following regularized increment term has been introduced, 

7^G?'^'+(l,2):=c2 / / e-^(l«il+l6l)|^^^2|^-"G+(6,6) (1-22) 

and 

nG'/'-{l,2):=cl [ [ e-^(l«^l+l«^l)|66|^'°G,-(6,6)- (1-23) 

The regularized Levy area satisfies the multiplicative property (ii) of the 
Introduction because the corresponding counterterm 

(5Z2'^(1,2) := ((57^G2'"+(l,2) - 6G^^''+ {1,2)) + {6TZG^''^- (1,2) - 6G^'''-{1,2)) 

(1.24) 

(given by an integral on the conical Fourier domain M? \ K^eg along the diag- 
onal 6 = ~6) is an increment (see Introduction). It satisfies the geometric 
property (iii) because (5Z-^'^(1, 2) is antisymmetric in 1 2, which follows in 
turn from the symmetry exhibited at the beginning of the proof of Lemma 
Ol 
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Theorem 1.1 For every a G (0, 1/2), TZB^fll, 2) is 2a -Holder uniformly 
in e. 

Proof. Let us first prove the Holder estimates (i) E\nBff {1,2)\^ < 
C\t — sl^"^. Similarly to the proof of point (ii) in Lemma 11.31 we apply 
Lemma [Tm (i) to Freg{u) = J^dW^UregiO^^^^ '^it^ (setting £, = £,i +^2) 

aregiO = -J' dT^^, " 6 1 ^""^^e-^d^^l+l^-^^ D , (1.25) 
K Jd 1?2 

where = e M | |6| > 1^-61, \C\ > Creg\C2\}. In particular, |C -^1 < 
^j^l^l, which implies |^2| < C'l'^l) so (assuming for instance 4 > 0) 

Vara,e5(0 < ^ T^^r'-'ie - 6r-'"ci6 

= ler'"^" / |u|l-2°|l-^z|-^-2" (i^Z = C'|er'-'" <00. 
Jl/2 

(1.26) 

For the rate of convergence, we rewrite e"^'''^^'^!^^!^ — e"''^!^^'^!^^'-' as the 
sum of two terms as in the proof of Lemma [1.3l (i). The second term may be 
studied using Lemma [L^ (ii) by considering the coefficient of dW^^(2); since 

^ < 1^2! < c]f^5 one gets the same estimates as for the proof of the Holder 
estimates (i), and hence a rate of convergence with exponent 4a. The first 
term involves the function F{e) = J (iW^^ajj(^i)e'"^ie~^l^il, with 

The estimates E|a^(^i)p < C|^|~^~^" is easy to get using the fact that 
Id + C2I > Creg\^2\ > CreglCil) which gives oncc again the same rate of 
convergence. 

□ 

Remark. Since the boundary terms Bj^^'^£l, 2)((9) are uniformly 2a~- 
Holder in e, one could simply have set TZG^f' (1,2)(9) = (which looks 
somewhat drastic). Our point of view is to try and regularize as little as 
possible. 
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2 Iterated integrals of higher order: a sketchy overview 



Combinatorics of Fourier normal ordering become non trivial starting from 
iterated integrals of third order; we shall concentrate on this case in this 
section, although some notions will be presented for the case of general 
iterated integrals. Detailed proofs should be found in 27j. 

Recall the two decompositions of Its{^i,^2) into 6Gfg{^i, ^2)+Itsi^i^ ^2){d) 
and (^G^(^2)Ci) + -^ts(?2) As already noted, the symmetry <-> ,^2) 
s ^ t maps to and I^{d) to I~{d), hence terms with indices it may 
be treated on an equal footing. This is no more the case for the integrals 
involved in B^^^. Namely, Fubini's theorem implies for instance (considering 
three among the six permutations of {1,2,3}, including the trivial one) 




Each of these may be represented as a finite sum of tree iterated integrals 
as we shall presently see. If T is a decorated rooted tree, i.e. a tree with 
a distinguished vertex vi called root, such that each vertex v G V{T) = 
{vertices of T} wears a label £{v) G {1, . . . , d}, and F = (F(l), . . . , T{d)) is a 
smooth d-dimensional path, then the integral of F along T is (denoting by 
v~ , V £ V{T) \ {vq} the unique ancestor of v, i.e. the unique vertex just 
below v) 

[h{r)]ts ■■= jyv^^^i^iv,)) jy dF.„^(^(t;2))... dT,^^^^^^^{l{v\y(r)\)) 

(2.4) 

where (ui, . . . , v\v{t)\) is any ordering of y(T) compatible with the tree par- 
tial ordering, i.e. such that < v for every v E y(T) \ {fi} (in other 
words, such that the indices of the vertices decrease while going down the 
branches towards the root). The definition extends easily to forests, i.e. to 
(finite) disjoint unions of trees (which are seen as a commutative product of 
the trees), by multiplying the tree iterated integrals corresponding to each 
connected component of the forest, and then (by taking linear combinations) 
to the algebra over M generated by decorated rooted trees, T, which 
is actually a Hopf algebra ^ ^ . 
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Now replace in eq. (I2.1l2.2l2.3p by - /" and /" by - Then 
eq. (j2.1l2.2l2.3|) may be represented resp. as Ij^ {B"^), Ij^ ^ {B^) — I12 2 {B^), 
/Ta,,(i?^)-/T3,2(i^^) (seeFig.EI). 



ST' 

i^^l iVl ij#3 iVi 

Figure 1: Example of iterated integrals. From left to right: Ti, T2.1, T2,2, Ta^i, T3^2- 



More generally, to each permutation o" of {1, 2, 3} corresponds a rewriting 
;ome finite sum, 

B?/(ii,i2,i3) =E5(^'jVTj(i?^), (2.5) 



of Bj^^(ii, 5 ^3) as some finite sum 



where g{cr,j) = ±1 is a sign. Note that each Ijj{B^) equals more explicitly 
(after permutation of the names of the ^-variables, — > ^o-(i)) 



^ s ^ s ^ s 

(2.6) 

where {vi,V2,V3) is the natural ordering of V(T'j) (compatible with the 
tree ordering) given by the position of the corresponding variable of inte- 
gration inside the iterated integral after applying Fubini's theorem, see eq. 
(j2.li2.2i2.3p or Figuredl and x - = t if is a root. Note also that the labels 

i 

on the trees are simply i{vj) = 

The idea of Fourier normal ordering consists in rewriting first B^]'^(ii, ^2, 13) 
as (letting S3 be the group of permutations of {1, 2, 3}) 

(^iJZJZai^^d) / / dW^c.i(vi))c^t^c„2(M2))^^^c.3(M3)) • 



5i|<l6l<l6 



V3 ) 

(2.7) 



13 



where {vi,V2,V3} are the vertices of TJ, so that innermost integrals bear 
highest Fourier indices. 

The next task is to get rid of divergences by adding some counterterms 
(or in other words, by discarding the contribution to the integral of adequate 
Fourier subdomains). If one wants the multiplicative property to remain 
true, this must be done in compatibility with lower-order counterterms (i.e. 
with the definition of the regularized Levy area at this stage). Hence one 
should have: 

(2.8) 

This identity has been generalized to tree integrals (see [H]). Namely, 
letting T be a tree and v range over all admissible cuts of T (we shall use 
the notation: v |= V(T)), i.e. over all non-empty subsets {vi, . . . ,vj} £ 
V{T) \ {0} (O=root of T) such that no pair {vi,Vj} C v is connected by 
going down or up the tree, then 

„T(-B^)]tu[/ii„T(-B^)]ns (2.9) 

v\=V(T) 

where R^T C T is the forest obtained as the union of branches lying 
above vi,...,vj (including the vertices vi, . . . ,vj), and L^T C T is the 
subtree obtained after removing these branches (see example on Fig. [2]). 
This identity called tree multiplicative property may be rephrased [28] 
by saying that [6I{B'^)]tus (viewed as a linear form on the Hopf algebra T) 
is the convolution of [I{B^)]tu and [I{B^)]us- 

^ • • ^ • • • 

1 13 2 1 12 3 

Figure 2: Admissible cuts of T2;2- v — {3}, {2} or {2,3} is also in this particular 
case the set of vertices of R^T. 

Now (see footnote after Lemma ll.Sp . one wants to set formally s = itioo 
to generalize the increment /boundary decomposition of the previous section. 
This results in the following definition of skeleton integrals: 



14 



Definition 2.1 (skeleton integrals) Let 

Jm. Jm. Jr 

e-^(iai+l61+l61) pi e^-^i^idx,, e^^^^^^dx,, p' e'^'^^^^dx,, 

(2.10) 

fsee eq. <fI76|) j. 

Formally [Sk/T(S")]t = [/T(-B^)]t,±ioo. As was the case for ^^'^'±(1,2), 
only the increment [5S]^If{B'^)\ts = [S]^If{B^)\t — [Sk/'f(i?^)]s makes sense. 
The tree multiplicative property yields then the tree skeleton decompo- 
sition 

[h{B')\tu = [SkJT(i?=)]t - [Skh{B')]u - [lL.T{B')\tu[^klR^T{B')\u. 

(2.11) 

This is an inductive formula, which yields Ij{B^) in terms of integrals 
or skeleton integrals of lower order. Once again, it may be interpreted as 
a convolution of characters of the Hopf algebra T, as the convolution of 
[Sk/(i?'^)]t with the inverse of [Sk/(i?^)]u (defined via the antipode) to be 
precise, see [28]. If e.g. T has 3 vertices, then each tree component of R^T 
has at most 2 vertices. A skeleton integral for a tree reduced to one vertex 
is simply 

/ dVFge-^l«l|e|5-" re'<du = [ dTyge-=l«le'"«^l^ (2.12) 

JR J Jr K 

(whose increments are those of B^ itself) and requires no regularization, 
whereas for a tree with two vertices, one gets 

/ dW^,{l) I dW^2(2)e-^(l«il+l«2|)|^^^2|^~° / e^'^-i^^dx^, [ ^'e'^^^^^dx^^ 
Jm Jr J J 

= [ dW^.il) [ dW^«,(2)e-^(l«^l+l«2l)|ei6l^-"G'^(ei,6) (2.13) 
Jr Jr 

(which is divergent in the limit e — > when a < 1/4, as proved in Lemma 
[Ll. 

The reader may easily check that formula (j2.1ip gives precisely the in- 
crement/boundary decomposition of the previous section when T has two 
vertices. 
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In order to get convergent quantities in compatibility with the regulariza- 
tion of second-order integrals, one should (i) regularize the new skeleton inte- 
grals of third order [SkIj{B'^)]t; (ii) replace [lL^T{B^)]tu and [SklR^j{B^)]u 
in the right-hand side of (j2.1ip by the corresponding regularized quantity of 
order (1 or) 2. 

In other words, regularization must be performed on each skele- 
ton integral of order > 2. 

Definition 2.2 (regularized skeleton integrals of order 3) Let, for T = 
Ti, ¥2^1, T2^2) TTa,! or ¥3^2 corresponding to one of the three above permuta- 
tions cr C S3, 

(2.14) 

where IRjFe„ C is the suhdomain of integration defined by: 







= {161 


< 


161 


< 


161; 


16 + 61,16 + 6 + 61 >c;,,|e3|} 


(^^) ] 




= {|6I 


< 


161 


< 


161; 


16 + 61 >c;,,J6|}; 


(Hi) ] 


M.reg 


= {|6I 


< 


161 


< 


161; 


16 + 6 + 61 >c^;eJ6|}; 


(iv) ] 




= {|6I 


< 


161 


< 


161; 


16 + 61 >c;,J6|}; 


(v) ] 




= {|6I 


< 


161 


< 


161; 


16 + 6 + 61 >c;,je3|} 



for some constant C^e„ G (0, 1). 



For the general definition of the cut Fourier domains of integrations, we 
refer the reader to [27j. Uniform Holderianity with respect to e has been 
defined before Lemma 11.31 

Theorem 2.1 Assume a < 1/3. Then the above regularized skeleton inte- 
grals are 3a~ -Holder uniformly in e. 

Proof (sketch). Let us just sketch the proof e.g. for [T^Sk/fi (-B^)]^, 
assuming i\ i2 is, say, ii = 1, 12 = 2, i^ = 3. The triple integral 
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J "i e'^^^i^Mxi,^ / "2 e''^'"2?2(ix^2 / e^^^s^^dxy^ writes (by straightforward 

computation) [1^3] [i (ga +6 )j[i^(^i^+g2 +g3 )] ' '^^^ ^"^^^ ^'^^ restriction of the 

domain of integration ensures that denominators are not too small. Apply 
Lemma 0(1) to F{u) = ciPVga(C)e'"? with (setting ^ = ^1 + ^2 + 6) 

on a domain _Df on which ^ < |^3| < Jr^, |^ — ^2 — Csl^^" < I'Csl^ ^ and 

^llJ < ct'^^r^-". Considering Var a(0 < § / ^^^d^sl^r'-'^", 
the integral over |^2| < iCsl contributes OdCsD) while the integral over ^ < 
161 < leads to gO(|e|i-^") = Oder'-'")- Hence E([,57^Sk/Tl(S^)]^s)' < 
C|t — s|^". By applying Lemma 11.41 (ii), a similar convergence rate may be 
proved. □ 

We may now finally define regularized iterated integrals of order 3 ac- 
cording to the above scheme. 

Definition 2.3 (regularized integrals of order 3) Let 

[nij{B')]tu{ii,i2,i3) = Yl 2^5(^T,j)[7^/TJ(i?')]^« (2.16) 

where g{(T,j) is as in eg. 112. 5\) . and by definition 

[TZhiB')]tu = [nskh{B')]t-[nskh{B')]u- mLMB')]tu[nskiii^j{B' 

v\=V{T) 

(2.17) 

The definition mimicks the previous unregularized tree multiplicative 
property for skeleton integrals, eq. ()2.1ip . All terms in it are convergent 
when £ ^ 0. There remains only to prove that the multiplicative/Chen (ii) 
and geometric/shuffle property (iii) of the Introduction are preserved by the 
regularization. In principle (ii) should be true because eq. (|2.17p is a mul- 
tiplicative property in itself; it must only be shown that this multiplicative 
property implies the original multiplicative property (ii) of the Introduction 
after summing over all permutations a, which is best done in a Hopf al- 
gebra language. The geometric/shuffle property is proved by showing that 
the regularized integration operator X^o- 5'(<^)i)^Sk/']r<T is a character of 
another Hopf algebra called shuffle algebra. Proofs should be found in [28j. 
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